%% Carbon Taxes Around the World: Cooperation, Strategic Interactions, and Spillovers
% IMF Economic Review
% Alessandro Moro and Valerio Nispi Landi
% Replication files
% This file avoids to run optimal_tax again to plot Figure 2 and 3, once you have sim_psi0, sim_psi05, sim_psi1




if psiX==0
load sim_psi0



ne=ne_psi0;
sp=sp_psi0;
coord=coord_psi0;
cequiv_sp=cequivz_sp_psi0;
cequiv=cequivz_psi0;
cequivz_sp=cequivz_sp_psi0;
cequivz=cequivz_psi0;
WW=WW_psi0;
WWz=WWz_psi0;
WWw=WWw_psi0;
DG= 0.00:0.0001:0.007;   
DGZ=0.00:0.0001:0.007; 
nn=length(DG);
nnz=length(DGZ);
DGend=0.02;
DGZend=0.02;
DG(nn+1:nn+1000*(DGend-DG(end)))=DG(end)+0.001:0.001:DGend;
DGZ(nnz+1:nnz+1000*(DGZend-DGZ(end)))=DGZ(end)+0.001:0.001:DGZend;
if LIN==1
DG=DG/10;
DGZ=DGZ/10;
end
if DAM==1
DG=10*DG;
DGZ=10*DGZ;
end

elseif psiX==0.5

load sim_psi05
ne=ne_psi05;
sp=sp_psi05;
coord=coord_psi05;
cequiv_sp=cequivz_sp_psi05;
cequiv=cequivz_psi05;
cequivz_sp=cequivz_sp_psi05;
cequivz=cequivz_psi05;
WW=WW_psi05;
WWz=WWz_psi05;
WWw=WWw_psi05;
DG= 0.00:0.0001:0.007;   
DGZ=0.00:0.0001:0.007; 
nn=length(DG);
nnz=length(DGZ);
DGend=0.02;
DGZend=0.02;
DG(nn+1:nn+1000*(DGend-DG(end)))=DG(end)+0.001:0.001:DGend;
DGZ(nnz+1:nnz+1000*(DGZend-DGZ(end)))=DGZ(end)+0.001:0.001:DGZend;
if LIN==1
DG=DG/10;
DGZ=DGZ/10;
end
if DAM==1
DG=10*DG;
DGZ=10*DGZ;
end

elseif psiX==1
load sim_psi1


ne=ne_psi1;
sp=sp_psi1;
coord=coord_psi1;
cequiv_sp=cequivz_sp_psi1;
cequiv=cequivz_psi1;
cequivz_sp=cequivz_sp_psi1;
cequivz=cequivz_psi1;
WW=WW_psi1;
WWz=WWz_psi1;
WWw=WWw_psi1;
DG= 0.00:0.0001:0.007;   
DGZ=0.00:0.0001:0.007; 
nn=length(DG);
nnz=length(DGZ);
DGend=0.02;
DGZend=0.02;
DG(nn+1:nn+1000*(DGend-DG(end)))=DG(end)+0.001:0.001:DGend;
DGZ(nnz+1:nnz+1000*(DGZend-DGZ(end)))=DGZ(end)+0.001:0.001:DGZend;
if LIN==1
DG=DG/10;
DGZ=DGZ/10;
end
if DAM==1
DG=10*DG;
DGZ=10*DGZ;
end

end